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Abstract 

We re-investigate the impact of isospin violation for extracting the s-wave irn scat- 
tering phase shift difference 6q(Mk) — 5%{Mk) from K — > tttt decays. Compared 
to our previous analysis in 2003, more precise experimental data and improved 
knowledge of low-energy constants are used. In addition, we employ a more ro- 
bust data-driven method to obtain the phase shift difference Sq(Mk) — ^{Mk) = 

(52.5 ± 0.8 exp ± 2.8 t hcor)°- 



1. If isospin were conserved the Fermi- Watson final-state interaction theorem would allow 
to extract the s-wave pion pion phase shift difference Sq(Mk) — 5 2 (Mx) directly from 
K — > 7T7T decay rates. However, the K — > tttt amplitudes are sensitive to isospin violation 
including electromagnetic corrections, especially the 1 = 2 amplitude A 2 . This is due to 
the large ratio A /A 2 ~ 22 (octet enhancement in nonleptonic weak decays) that enhances 
isospin- violating corrections to A 2 . 

It has been a long-standing problem to reconcile the phase shift difference extracted 
from K — > 7T7T decays with other determinations of pion pion phase shifts. This problem has 
become especially acute after the precise determination of tttt phase shifts from combining 
dispersion theory with chiral perturbation theory [I]. Our previous analysis of isospin 
violation in K — > tttt decays [2] (similar results were obtained in Ref. [3]) led to S (M K ) — 
S2(M K ) = (60.8 ± 2.2 cxp ± 3.1 t heor)°, substantially bigger than the dispersion theoretical 
result [1] 5 (M K ) - S 2 (M K ) = (47.7 ± 1. 5)°. 

We have decided to reanalyse the problem for several reasons. 

• The experimental situation has substantially improved since 2003 for both the K + 
and K s lifetimes and for the branching ratios of K — > nn decays [4-6]. It was already 
observed in Ref. [7] that the new experimental information reduces the phase shift 
difference by more than three degrees (with higher statistical significance), bringing 
it closer to the dispersion theoretical value. 

• New information has also become available on some of the low-energy constants 
(LECs) involved, both in the strong [8] and in the electromagnetic sector [9]. The 
effect on the phase shift difference is difficult to quantify but it is definitely smaller 
than the uncertainty assigned to the LECs in Ref. [2]. As one example, the new 
estimates of electromagnetic LECs [9] lead to a shift of the lowest-order it — rj 
mixing parameter from 1.06- 10~ 2 to 1.29- 10~ 2 [10]. This shift reduces the phase 
shift difference by 0.2°. 

• As will be detailed below, the theoretical analysis of Ref. [2] can also be improved by 
relying to a lesser extent on the NLO calculation of K — > tttt decay amplitudes. Using 
less information potentially increases the uncertainty but this will be compensated 
by a less biased comparison with the data. The resulting estimate of isospin violation 
is expected to be more robust than the original estimate and it leads to a further 
decrease of the phase shift difference by nearly two degrees. 

• Based on a method proposed in Ref. [11] making use of the optical theorem, effects of 
0{e 2 p A ) were partially accounted for in Ref. [2]. Since a complete calculation of such 
effects is beyond present capabilities, we have decided to keep track of the associated 
theoretical uncertainty, but without including the partial corrections in the mean 
value for the phase shift difference. 
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2. We start by recalling the procedure of Ref. [2] for extracting the phase shift difference. 
The amplitudes A + _ = A{K° -> tt+tt-), A 00 = A(K° -> 7r°7r°) and A +0 = A(K + -> tt+tt ) 
are parametrized as 

v2 

A 00 = A e <xo - v^A 2 e ixa (1) 

o 

^+0 — 2 2 

In the absence of CP violation, the amplitudes A , A 2 , A 2 are real and positive by definition. 
In the isospin limit, A 2 = A^ and the phases xi coincide with the strong rnr phase shifts 
5i at the kaon mass. 

To NLO in the chiral expansion, the phases xo, X2 cannot be calculated reliably: the 
resulting phases are substantially too small. To obtain the strong phase shift difference 
5 (M K ) — 5 2 (M K ), a two-step procedure was employed [2]. Using the NLO expressions for 
the absolute values A , A 2 , A£, the lowest-order couplings g 8 , g 2 7 were determined from the 
experimental rates together with the phase shift difference Xo ~ X2- Comparing the NLO 
amplitudes with and without including isospin violation, the differences 

V = Xi-Si(M K ) (7 = 0,2) (2) 

were calculated to obtain the final phase shift difference 8q{Mk) — S 2 (Mk) in the isospin 
limit. 

There are two related potential pitfalls associated with this procedure. Although the 
theoretical NLO expressions for the phases xi cannot be trusted they enter the dispersive 
and absorptive parts of the amplitudes A , A 2 , A\ implicitly. Moreover, although both the 
Xi and the 5j cannot be calculated reliably to NLO in the chiral expansion the differences 
7/ were assumed to be trustworthy. Both steps are therefore subject to a theoretical bias 
that is difficult to control at the level of accuracy considered. 

The main idea of the alternative procedure proposed here is to use only the isospin- 
violating parts of the NLO amplitudes as theory input and to determine So(Mk) — 5 2 (Mk) 
directly from the data. In contrast to the chiral corrections for the full amplitudes, the 
isospin-violating corrections are much smaller and therefore less subject to the bias dis- 
cussed in the previous paragraph. 

The amplitudes are now parametrized as 

A + _ = A e iSo ^ + ^-A 2 e iS ^ M ^ + AA™_ 

v2 

A 00 = A e i5 °( M ^ - V2A 2 e i52 ^ + AA^ (3) 
A +0 = ^A 2 e iS ^ + AA%. 

All isospin violation is contained in the amplitudes AA 1 ^, AA 1 ^, AA 1 ^. They can be 
extracted from the NLO amplitudes of Ref. [2]. Since isospin violation was neglected in 
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the 27-plet amplitudes because of the AI = 1/2 rule [2,3] the amplitudes AA™ (n = 
H — , 00, +0) scale linearly with the lowest-order octet coupling g s . They also depend on 
higher-order LECs in addition to the loop contributions. 

The moduli of the amplitudes in the isospin limit0 are denoted as A , A 2 . We will not 
use the theoretical expressions for these amplitudes but instead determine them together 
with the phase shift difference S (M K ) — 5 2 (M K ) directly from the rates. For this purpose, 
we write the moduli of the amplitudes as 



\A 00 \ 
l^+o I 



]4 _|_ J_~^ 2e i{&2{M K )-& {M K )) _|_ ^IB e -iS (M K ) 

A - V2A 2 e ^{M K )-8 (M K )) + aa ib e -i6 (M K 



-Ao ^(MMKysoiMK)) , AA m 
2 



+o' 



00 

i<5 (M x ) 



(4) 



In order to determine A , A 2 and Sq(Mk) — 5 2 {Mk) from the three rates, we therefore also 
need the 1 = phase 5 (M K ) as input. We use the value obtained in Ref. [1]: 



6 (M K ) = (39.2 ±1.5)° . 



(5) 



From the structure of the moduli (j4j) one may already anticipate that the precise value of 
5q(Mk) will have little impact on the phase shift difference. 



3. We use the same experimental input as in Ref. [7], which is reproduced in Table [U 



Parameter 


Value 


Correlation 


Reference 


BR(K S -> 7T + 7T-(7)) 

BR(K S -> 7r°7r°) 


0.69196(51) 
0.30687(51) 


-0.9996 


[4] 




0.08958(5) ns 




[5] 


BR(K + -> 7i + 7i°) 


0.2064(8) 
12.384(19) ns 


-0.032 


[6] 



Table 1: Experimental input taken from Ref. [7]. 

In addition to the experimental input, we need the isospin- violating amplitudes AA^ B . 
With the central values of the various LECs and displaying explicitly the linear dependence 
on the octet coupling g s , we find 

AA™_ = g 8 [2.25 - 0.83 «] • KT 10 GeV 

AA™ = #8 [-0.58- 2.77 1] • 10~ 10 GeV (6) 
AA% = g 8 [-2.12 - 1.10 i] ■ 10" 10 GeV . 

J Wc define the isospin limit in terms of the neutral meson masses [2]. 
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Updating the fit in Ref. [2] with the new data, we get a mean value g$ = 3.6 (see also 
Ref. [7]). We are going to assign a 20% uncertainty to g$, much bigger than the fit error. 

Using the isospin-violating amplitudes ([HD with g 8 = 3.6 and S (M K ) = 39.2° [1], the 
experimental rates in Table [1] give rise to 

A = 2.7030(8) • 10" 7 GeV 

~A 2 = 0.1249(3) • 10" 7 GeV (7) 
S (M K ) - 5 2 {M K ) = (52.54 ± 0.83)° . 

The errors are purely experimental but they take the correlations in Tabled] into account. 
The resulting correlations for the fitted quantities are small, at most —16% for the corre- 
lation between A 2 and S (M K ) — 52(M K ). 

The octet enhancement in K — ► -kit decays is characterized by the amplitude ratio 

= = 2L63 ( 4 ) > ( 8 ) 

again with experimental error only. 



4. For assessing the theoretical uncertainties, we concentrate on the phase shift difference. 

The uncertainty induced by the error of Sq(Mk) in (jHJ) can be disposed of quickly. 
Varying the 1 = phase shift as 6q(Mk) = (39.2 ± 3.0)° (2 a) affects the phase shift 
difference only in the second decimal place for 5 (M K ) — S 2 (M K ). Of course, this has to do 
with the smallness of the isospin-violating amplitudes ([6]). The uncertainty in the phase 
shift difference is therefore completely negligible. 

The overall scale of the isospin-violating amplitudes (jSJ) is determined by the lowest- 
order coupling g 8 . We assign a generous error of 20%, i.e., g 8 = 3.6 ± 0.8. As already 
emphasized, this uncertainty is much larger than the fit error [2, 7] but it is meant to 
account also at least partially for effects of 0(e 2 p 4 ) (see below). Varying g s by 20% gives 
rise to an uncertainty ±1.1° for the phase difference. 

As in Ref. [2] , we estimate the uncertainty associated with the various LECs by varying 
both the short- distance scale for the Wilson coefficients (0.77 GeV < /xsd < 1-3 GeV) 
and the chiral renormalization scale (0.6 GeV < u x < 1 GeV). The central values in (J7|) 
correspond to //sd — 1 GeV, v x = 0.77 GeV. As already mentioned, the changes in the 
LECs from 2003 till today are well within the range expected from varying the two scales. 
It turns out that the error associated with the short-distance scale is asymmetrical: the 
phase difference happens to be minimal for /xsd — 1 GeV and the increase from varying 
/xsd is at most 0.5°. From the dependence on the chiral renormalization scale, the error 
for the phase difference is bigger and nearly symmetrical around the central value for v x = 
0.77 GeV: the phase difference varies by ±1.2°. 

Finally, we consider higher-order effects of 0(e 2 p 4 ). In Ref. [2] we made use of the optical 
theorem and tttt scattering amplitudes to 0(e 2 p 2 ) [12,13] to estimate the absorptive part of 
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the K — > mi amplitudes to 0(e 2 p 4 ). Based on that we quantified the 0(e 2 p 4 ) corrections 
to 72 (and to 5 (M K ) — 5 2 {M K )) to be +2.6° and included them in our final estimate. 
However, since a complete calculation of K — > txtx to 0(e 2 p 4 ) is not available, here we 
assume a more prudent attitude of not including this correction in the central value. We 
rather take this result as a first measure of the size of higher-order corrections. We have 
also reached a similar conclusion on the size of these effects with no reference to the optical 
theorem analysis. Within the method introduced in this letter, we have parametrized 
higher-order corrections via three different scale factors multiplying the isospin-breaking 
amplitudes AA^_ 00 +0 . We have repeated the fit for different choices of the scale factors 
ranging independently between 0.5 and 1.5, finding that the output 8q{Mk) — S 2 (Mk) 
changes by at most 2.2°, which will be used as our final estimate of higher-order effects. 

Altogether, our new procedure for confronting theory with experiment leads to the 
following final value for the phase shift difference in the isospin limit: 

5 (M K ) - 5 2 (M K ) = (52.5±0.8 cxp ±l.l 98 ^ SD ± 1.2 X ± 2.2 0(eV) )° 

= (52.5 ± 0.8 cxp ± 2.8 thcor )° . (9) 



5. If isospin violation is neglected (except in the physical pseudoscalar masses for phase 
space), i.e. taking AA 1 ® = 0, the fit to the experimental rates gives 



[6 (M K ) - 8 2 {M K )} 



Isospin 



(47.3 ± 1.0) c 



and 



A 

A, 



22.41(5) . 



(10) 



Isospin 



The substantial experimental improvements achieved with the most recent kaon data have 
reduced the phase shift difference from the value (48.6 ± 2.6)° obtained in 2003 in the 
isospin limit [2] (and about 58° some 30 years ago [5]). The value (flUl) would be in perfect 
agreement with determinations from tttt scattering data: 



[S (M K ) - 5 2 {M K )l 



(47.7 ± 1.5)° 
(50.9 ± 1.2)° 



[1] 
[14] 



(12) 



However, owing to the large A /A 2 ratio in K — > 2ti decays, isospin-breaking corrections 
to the dominant AI =1/2 amplitude generate sizeable contributions to A 2 [compare the 
results (jHJ) and (TTTj) ] . modifying also the amplitude phases. Our 2003 analysis of isospin 
breaking in K — > 2tt decays concluded that these effects increase the phase shift difference 
significantly. Including some 0(e 2 p A ) corrections through the optical theorem, we found 
the result 5 (M K ) - 5 2 (M K ) = (60.8 ± 2.2 exp ± 3.1 t heor)° [2]. The large difference with the 
7T7T determinations ffl2|) has been a pending puzzle since then. 

In this letter we have reanalysed the K — > 2tt phase shift determination, taking advan- 
tage of the improved experimental situation. Moreover, we have modified the theoretical 
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analysis in order to be less sensitive to theoretical uncertainties. Since the absorptive con- 
tributions that generate the strong amplitude phases start to appear at the one-loop level, 
a NLO theoretical calculation of the amplitudes only provides the phase shifts at leading 
order, which are therefore subject to large uncertainties. To minimize theoretical errors, 
we have only used as theory input the calculation of the isospin-breaking contributions 
AA 1 ®. In this way, we can determine all other quantities (the phase shift difference and 
the amplitudes in the isospin limit) directly from a fit to the data. The residual theoretical 
uncertainties associated with the isospin-breaking contributions have been estimated con- 
servatively in two different ways, as explained in the previous section. Our final result in 
Eq. ([9]) is lower than our 2003 determination. Both the new data (as already observed in 
Ref. [7]) and the different treatment of theory input tend to lower the resulting value: 3.3° 
from experiment and altogether 5° from theory. This updated determination of the phase 
shift difference from K — > 2n decays turns out to be in agreement with the tttt results in 
f ll2p . although with a larger uncertainty. 
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